Evidence of quantum phase slip effect in titanium nanowires 
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Electron transport properties of titanium nanowires were experimentally studied. Below the effec- 
tive diameter < 50 nm all samples demonstrated a pronounced broadening of the R(T) dependencies, 
which cannot be accounted for thermal flcutuations. An extensive microscopic and elemental anal- 
ysis indicates the absence of structural or/and geometrical imperfection capable to broaden the the 
R(T) transition to such an extent. We associate the effect with quantum fiucutuations of the order 
parameter. 

PACS numbers: 74.25. F-, 74.78.-w 



Since the early years of experimental studies in su- 
perconductivity it has been noticed that the supercond- 
cunting transition R(T) has always a finite width. Very 
often the broadening can be accounted for sample in- 
homogeneity. However, soon it became clear that, at 
least in low dimensional samples, the transition width 
remains finite even with the refined material purity and 
improved fabrication. The effect has been attributed to 
fluctuations typically more pronounced in objects with 
reduced dimensionality. The finite resistance R(T) ~ 
exp (— Fo/ksT) at a temperature T below the critical 
temperature T c of a quasi-one-dimensional superconduct- 
ing channel with cross section a has been explained by 
the thermal fluctuations of the order papramctcr: the 
so called thermal activation of phase slips (TAPS)fi^. 
Here the condensation energy Fo ~ B^a of the small- 
est statistically independent volume £er, where £ is the 
supercondcuting coherence length and B c is the criti- 
cal magnetic field, competes with the thermal energy 
ksT. The effect manifests itself only sufficiently close 
to the critical temperature, and in extreemly homoge- 
neous samples with micrometer-size diameter (e.g. pure 
whiskers) leads to the experimentally observable width 
of the R(T) transition of about few mK— In less 
homogeneous objects (e.g. lithographically fabricated 
nanowires) separation of the impact of the thermal fluc- 
tuations from the trivial inhomogeneity-determined R(T) 
broadening is rather problematic 6 . Nevertheless with de- 
velopment of nanotechnologyi it became clear that in 
extreemly narrow superconducting wires, with diame- 
ters ~10 nm, the shape of the R(T) transition by no 
means can be explained by sample inhomogeneity or/and 
thermal fluctuations^. The effect has been attributed 
to quantum fluctuations, also called - quantum phase 
slips (QPS) - and has been observed in a rather lim- 
ited number of experiments studying the transport prop- 
erties of ultra-narrow nanowires made of various super- 
conducting materials: amorphous MoGe~r^r^; In and 
In - P&i^jii; Being not explicitly at- 

tributed to quantum fluctuations, the anomalous R(T) 
broadening has been reported in Phl^ and amorphous 
InO;J^. Though the subject of a reliable experimental 



FIG. 1: (Color online). SPM image of typical part of a 
titanium nanowire. Horisonal plane indicates the interface 
between the metal and the sputtered Si substrate. The ex- 
periemental error in definition of the interface position pro- 
vides the main contribution to uncertainty in determination 
the sample cross section. Inset: profile of the top part of the 
sample. 



confirmation of the quantum fluctuation phenomenon in 
quasi-one-dimensional superconductors is still under de- 
bates, there is a consensus in the scientific community, 
that if exists, it should be observed in extreemly narrow 
samples with characteristic diameters ~10 nm. Unfortu- 
nately, at these scales an independent and reliable analy- 
sis of a nanostructure homogeneity is rather problematic. 
The uncertainty leaves a room for a critically-oriented sci- 
entist to attribute the experimentally observed deviations 
of the R(T) superconducting transition shape from the 
well-established TAPS model to sample inhomogenities: 
structural (impurities, grain boundaries) or geometrical 
(constrictions). In this paper we demonstrate that with 
proper selection of material (superconducting titanium) 
the obvious deviations of the R(T) transition shape from 
the TAPS model become pronounced already at scales 
< 50 nm, which drammatically simplifies both the sam- 
ple fabrication and makes the structural analysis more 
reliable compared to the 10 nm case. 

So far a rather limited number of models has been 
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proposed to describe the impact of quantum fluctua- 
tions on transport properties of quasi-one dimensional 
superconductors^. Followin g 20 , 21 , 22 one may describe 
the QPS contribution to the effective resistance of a su- 
perconducting wire with length L and cross section a as: 

R(T) * b ^ PS cxp(-2S QPS ) (1) 

where b is an unimportant constant which remains 
the same for all samples; A(T) and £(T) are the 
temperature-dependent superconducting energy gap and 
coherence length, respectively. The QPS action Sqps = 
A(R Q /R N )(L/£), where R Q = h/2e = 6.47 kfl is the 
'supercondcuting' quantum resistance and Rn is the wire 
resistance in the normal state. Constant A is of the order 
of unit and, unfortunately, cannot be determined more 
precise within the model^S^^ 2 -. The mean free path 
I and the constant A are the two fitting parameter of 
the QPS model. Note that the mean free path is not a 
truly free parameter and with acceptable accuracy can be 
estimated from the normal state resistivity as the prod- 
uct IpN is the material constant. Two other parame- 
ters of the model - the critical temperature T c and the 
normal state resistance Rn - are deduced from the the 
experimental R(T) dependencies. Recovering that for a 
'dirty limit' superconductor I <C £ the coherence length 
£ ~ \/I^o, where the 'clean' coherence length £o ~ hvp/A 
and vp is the Fermi velocity, one may conclude that with 

1/2 

the exponential accuracy R(T) ~ exp(— aT c <t/pn)- As 
the Fermi velocity does not vary much between con- 
ventional supercondcutors being of the same order as 
vp ~ 1.79 X 10 6 m/s for titanium, the constant a should 
be basically material-independent. So far in the rather 
limited number of experiments claiming the observation 
of QPS phenomen a 9 1 the efforts were mainly concen- 
trated on reducing the sample cross section a. However, 
the material issue was largely ignored: for the smallest 
obtainable dimension a, obviously limited by available 
fabication capabilities, the ultra-low temperature super- 
condcutors with high normal state resistivity are of ad- 
vantage. 

In this paper titanium has been selected as the suit- 
able material for demonstration of the QPS phenomenon. 
The crtitical temperature T c of the titanium nanowircs 
is below 400 mK. The low-temperature resistivity Pn 
is significantly larger than for the majority of single- 
element supercondcutors and varies from ~ 1.0 x 10~ 6 
ttxm for the 2D films to ~ 3.2 x 10~ 6 ttxm for the 
sub-30 nm nanowires. Utilizing for titanium the product 
ipN — 10 x 10~ 16 Slxm 2 , which slightly varies from dif- 
ferent literature sources^, one gets that in our nanowires 
the mean free path (. is of the order of 1 nm, which cor- 
relates with the independent transmission electron mi- 
croscope (TEM) analysis indicating the corresponding 
size of the defect-free areas. From the technological 
point of view titanium is an easy-to-deal material. The 
nanowires with different lengths L between 1 /jm and 100 
/im were fabricated using conventional lift-off technique: 
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FIG. 2: (Color online) Resistance vs. temperature for the 
same titanium nanowires with length L= 20 /im and pro- 
gressively reduced effective diameter %fo indicated in the plot 
and specified with accuracy ±2 nm. Inset: normal state re- 
sistance per unitlength vs. inverse cross sections for several 
sample swith different length L indicated in /^m. 



soft-mask e-beam lithography followed by e-gun evapo- 
ration at a residual pressure ~ 10 _9 mBar on naturally 
oxidized Si/ SiO x substrate. The mentioned difference in 
normal state resitivity pjv between the wide films and the 
thin nanowires presumably originates from the pressure 
gradient between the bottom of the narrow groove in the 
resist mask and the rest of the vacuum chamber. In the 
former case, being an effective getter material, Ti 'ab- 
sorbs' the residual gas and the not-completely-evacuated 
organics leading to formation of a dirtier sample. Af- 
ter the analysis with the scanning electron and scanning 
probe microscopes (SEM and SPM, respectively), the 
samples showing no obvious defects were cooled down in 
3 He 4 He dilution refrigerator down to temperatures T ~ 
50 mK. Conventional four-probe DC and low-frequency 
(< 20 Hz) lock-in AC techniques were used to measure 
the R(T) dependencies. Special care was taken not to 
overheat the samples. Very low excitation currents down 
to 30 pA were used to ensure the linear response. The 
refrigerator and the front-end battery powered analogue 
pre-amplifiers were located inside the electromagnetically 
shielded room being connected to the external electronics 
via low-pass RC filters. An extensive multi-stage RLC fil- 
tering of the signal lines inside the refrigerator effectively 
reduced the impact of the noisy electromagnetic environ- 
ment. In a separate experiment using the same measur- 
ing set-up, the increase of the effective electron temper- 
ature T e deduced from the shape of the I — V depen- 
dence of a normal-insulator-supercondcutor junction 24 , 
was found to be at the level ST e < 20 mK above the base 
temperature T— . 

Both the co-deposited 2D films with characteristic 
thickness d ~40 nm and the as-fabricated relatively 
'wide' w > 60 nm nanowires showed an abrupt supercon- 
ducting transition at a critical temperature T c ~ 400 mK 
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FIG. 3: (Color online). Normalized resistance vs. temper- 
ature for three different nanowires of length L. The effective 
diameters yfa are indicated in the plot in nm and are specified 
with accuracy ±2 nm. 

(Figs. 2 and 3). The experimentally measurable width of 
the R(T) transition 5T C ~ 20 mK of these 'thick' samples 
can be qualitatively understood in terms of the TAPS 
model and the inevitable mild inhomogeneity, e.g. few 
percent variation of the cross section 6 . After each session 
of the R(T) measurments the samples were subjected to 
the low-energy ion etching to reduce the wire effective 
diameter y/a = \fd~w with steps as small as ~1 n m 26 , 27 . 
Due to the re-deposition of the sputtered material the in- 
evitable imperfections are smoothed-out resulting in the 
surface roughness of the processed samples approaching 
± 1 nm (Fig. 1). The method enables the study of 
a size-dependent phenomenon on a same structure with 
progressively reduced characteristic dimension elimintat- 
ing the artefacts of the samples fabricated in different ex- 
perimental runs. Further decrease of the nanowires cross 
sections leads to (i) reduction of the critical temperature 
T c> and (ii) broadening of the R(T) transition (Figs. 2 
and 3). The first effect is well-known for low dimensional 
superconductors^, though the origin of the phenomenon 
is still under debates 2 ^. As the size dependence of T c is 
observed not only in nanowires, but also in wide 2D films, 
we believe that the effect is not related to the essentially 
ID phenomenon under discussion - phase slips (thermal 
or quantum). 

Homogeneity of samples is the key point in interpreta- 
tion of experimental R(T) data within any model derived 
for a homogeneous superconducting channel of uniform 
cross section a. Critically oriented reader might always 
argue that the observed broad R(T) transition is the re- 
sult of a trivial sample inhomogeneity: either structural 
(e.g. local variation of the critical temperature along the 
wire) or geometrical (e.g. constrictions). The two cases 
should be analysed separately. Let us first consider the 
structural inhomogeneity. Obviously all real samples do 
have a certain level of structural inhomogeneity originat- 
ing from various sources: grain boundaries, finite size of 



FIG. 4: (Color online). V(J) characteristics of the three 
neighbouring parts with equal lengths L — 20 /im of the same 
nanowire with the effective diameter \fd = 38 ± 2 nm . Above 
the critical temperature T c ~ 300 mK within the scale of the 
image the Ohmic V(I) characteristics (o) are quantitatively 
indistinguishable between the different parts. Below the crit- 
ical temperature the V(I) dependencies are also very similar 
(■,»,A). Note the absence of the true zero resistance state 
below the 'critical current' value ~ 1.2 nA. Inset shows the 
layout of the sample and the measurement. 



the sample, proximity effect at the interface with wider 
parts of the structure, etc. We would like to note that 
all our samples, including the thinnest wires with a typi- 
cal sheet resistance i?n < 200 fl, are still comfortably on 
the metal side of a 'dirty' titanium. Formation of a net- 
work of weakly coupled metal grains and the correspond- 
ing Coulomb effects have been observed in deliberately 
oxidized titanium films with the sheet resistance Rq ex- 
ceeding few k O 29 ; 30 . All our samples above the critical 
temperature demonstrate I — V dependencies without 
any non-linearities, which otherwise might indicate the 
existence of weak links (Fig. 4), and the dependence 
of the normal state resistance Rn on diameter follows 
the expected Ohm's law (Fig. 2, inset). If the trans- 
port properties of the nanowires would be determined by 
the presence of the weak links, then at low temperatures 
the non-linear I — V characteristics should be sample- 
dependent. For example, various parts of a long multi- 
terminal nanowire should demonstrate some 'fingerprint' 
features related to formation of weak links particular for 
each sub-sample, which is not the case (Fig. 4). 

An extensive TEM analysis (Fig. 5) cannot reveal any 
suspicious structural imperfections inside the metal ma- 
trix: the material bulk looks exactly the same for the 
sputtered and for the non-sputtered samples. The poly- 
crystalline nanostructures consist of compactly packed 
grains with the average size of defect-free area ~ 3 nm, 
which correlate well with the the best-fit value I ~ 1 nm 
used in calculations. The elemental depth profiles (Fig. 
6) were determined by means of Time-of-Flight Elastic 
Recoil Detection Analysis (TOF-ERDA) 3 ! using 8.015 
MeV 35 C1 4+ incident ions. The areal density of oxygen 
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at the surface (1.7xl0 16 atoms/cm 2 ) and at the inter- 
face (5.5xl0 15 atoms/cm 2 ) corresponds to the thickness 
of about 1.9 nm and 0.8 nm, respectively, for the sur- 
face Ti02 and the boundary with Si02/Si with the cor- 
rersponding density 4.0 g/cm 3 and 2.2 g/cm 3 . The thick- 
ness of the material with the high concentration of oxy- 
gen correlates well with the high-resolution TEM anaysis 
data (Fig. 5). The bulk concentration of oxygen ~ 0.4 at. 
% inside the titanium matrix was determined by compar- 
ison of the experimental energy spectra with the ones ob- 
tained by Monte Carlo simulations 3 ^. Given that the av- 
erage microcrystall size is ~ 3 nm, one can easily estimate 
that inside the titanium bulk there is less that one oxygen 
atom per defect boundary. Concentration of other than 
oxygen elements inside the titanium matrix was found to 
be even smaller. The observation eliminates the possi- 
bility of the weak link(s) formation due to non-metallic 
grain interfaces capable to block the metal-to-metal su- 
percurrent. To summarize, all available at our disposal 
methods of analysis - SEM, TEM, SPM and TOF-ERDA 
- give us a confidence to state that (i) our nanostructures 
are as homogeneous, as a conventional thin film titanium 
can be; and (ii) - what is even more important - our 
method of reduction of the nanowire cross section by the 
low energy ion sputtering does not introduce new defects. 
At acceleration energies ~ 1 keV the penetration depth 
of the Ar + ions inside the titanium matrix is below 3 nm 
making the method virtually non-destructive: the thick- 
ness of the ion-damaged layer is comparable to the thick- 
ness of the naturally grown oxide. Hence, if there are 
some inevitable 'intrinsic' structural defects, they cannot 
appear in thinner samples contributing to broadening of 
the R(T) transition. If to assume that there exist some 
mysterious and undetectable mechanism which degrades 
the wire homogeneity with reduction of its diameter (e.g. 
recovery of hidden in the bulk cavern(s) or/and grain 
boundaries, blocking the supercurrent), it is reasonable 
to assume that the mechanism should manifest itself in- 
dividually for each particular sample leading to a unique 
'fingerprint' on the R(T) and the I — V dependencies. 
However, a large statistics of samples states the oppo- 
site: the shape of the broadened R(T) transitions is uni- 
versal and reproducible for different samples of the same 
effective diameter (Fig. 3). 

Now let us turn to the alternative critics dealing with 
geometrical inhomogeneity: the inevitable variation of 
a nanowire cross section. One might argue that the 
R(T) broadening originates from the well-known TAPS 
mechanism: due to the strong exponential dependence 
R(T) ~ exp(— a) the contributions coming from the parts 
with different cross sections ai oca i wash-out the otherwise 
sharp R(T) dependence. First, we would like to note that 
with reduction of the wire dimensions (by ion sputtering) 
the shape of the R(T) transition becomes more consis- 
tent between the samples with close values of the cross 
sections (Fig. 3). The observation supports the earlier 
statement that the ion beam etching polishes the surface 
making the sample geometrically more uniform: the as- 
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FIG. 5: (a) Bright field low-resolution TEM image of the 
cross section of a typical 35 nm thick Ti film on Si substrate, 
(b) Bright field high-resolution TEM image of the Ti / Si 
interface. Inset: Fast Fourrier Transform (FFT) indicates the 
single crystall nature of the metal grains forming the film. 



fabricated 'wide' nanowires have larger variation of the 
cross section a compared to the thinner ('polished') sam- 
ples. Here we would like to stress, that the uncertainty of 
the effective diameter specified in Figs. 2 and 3 reflects 
the experimental error in determination of the cross sec- 
tion and not the actual roughness of the surface, which 
for the multiply sputtered samples does not exceed ±1 
nm (inset in Fig.l). The error mainly comes from the un- 
certainty in determination of the position of the interface 
between the metal and the sputtered substrate (Fig. 1). 
In SEM the contrast between the two light materilas - Ti 
and Si - is not sufficent to determine the postion of the 
interface with accuracy better than ±2 nm. While the 
tip deconvolution effect, typical for SPM analysis of cs- 
sentailly 3D nanostructures, leads to the basically same 
uncertainty. Summarizing, we would like to state that 
by standards of the modern nanotechnology and micro- 
scopic analysis our nanowires are so 'large' that it is al- 
most impossible to overlook a pronounced constiction. 
Even in the worst case scenario, in the sputtered samples 
the variation the cross sections along the wire is below 
±10%. The contribution of such a moderate geometrical 
imperfection on the shape of the R(T) transition deter- 
mined by thermal fluctuations has been analysed 6 : by 
no means it can account for the experimentally observed 
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FIG. 6: (Color online). Experimental and Monte Carlo simu- 
lated energy spectra for oxygen recoils from the 33 nm thick 
Ti film and Si substrate. In the simulation the surface and 
interface oxygen contents were fixed, while the concentration 
in the film was varied. The best fit can be obtained with 0.4 
at.% concentration. The full depth profile of oxygen is shown 
in the insert. 



pronounced broadening of the R(T) dependencies. One 
may come to the same conclusion just analysing Fig. 2: 
the reduction of the wire average diameter from 58 nm 
to 33 nm leads to a negligible broadeneing of the TAPS 
fits (dashed lines) compared to the experimental data 
(symbols). For sufficiently narrow nanowires no realistic 
set of fitting parameters of the TAPS model can explain 
the broad experimental R(T) dependencies (e.g. Fig. 2, 
symbols). 

On the contrary, the QPS model provides reasonable 
agreement with the experiment (Fig. 2, solid lines). In 
simulations of the theoretical QPS-governed R(T) depen- 
dencies (Eq. [T]) for each sample (cross section) the best 
fit critical temperature T c corresponds to the onset of su- 
perconductivity, the mean free path £ ~ 1 nm and the 
numerical constant A ~ 0.3 were kept as free parameters. 
At temperatures T <^T C the negative magnetoresistance 
of about few percent has been observed in the thinnest 
samples at very small magentic fields < 3 mT. Similar ef- 
fect has been earlier reported in ultra-narrow lead^ 3 . and 
aluminumii nanowires. The origin of the pheonomenon 
is not yet clear. One alternative 3 ^ employs possible for- 
mation of a charge imbalance region accompanying each 
phase slip event. This non-equilibrium region, if exists, 
would provide dissipation outside the core of a phase slip. 
Within a certain range of (small) magnetic fields the cor- 
responding Ohmic contribution can be suppressed by the 
magnetic field more effective than the superconducting 
gap, resulting in the negative magnetoresistance. How- 
ever, so far the validity of the charge imbalance concept 
has been only demonstrated at temperatures sufficiently 
close to T c and its applicability to QPS is by no means ob- 
vious. The charge imbalance scenario responsible for the 
negative magnetoresistance 3 ^ still requires a solid theo- 
retical justification. 

It should be noted that the both models (QPS and 
TAPS) were derived assuming that the phase slips are 



the 'rare' events. In other words, the justified comparison 
with experiment is valid only in the limit R(T) <C i?j\r. 
The QPS effect is an essentially low temperature phe- 
nomenon when quantum fluctuations of the order param- 
eter dominate over the thermal fluctuations, with the lat- 
ter being important only sufficiently close to the critical 
temperature T — > T c . Contrary to the QPS, extrap- 
olation of the TAPS mechanism down to lower temper- 
atures violates the Ginzburg criterion (T c — T) /T c <Cl 
of the model applicability 3 ^. This is an additional ('the- 
oretical') argument why the broad experimental R(T) 
transitions (Figs. 2 and 3) cannot be explained in terms 
of thermal fluctuations^, even if one would assume in 
those samples the presence of extended and unrealisti- 
cally narrow constrictions y/a ~1 nm (overlooked in all 
microscopes!). 

In conclusion, we have studied titanium nanowires 
with progressively reduced cross sections. An extensive 
microscopy and elemental analysis revealed no obvious 
structural or geometrical imperfections. Neither the nor- 
mal state, nor the supercondcuting transport properties 
provided any signature of a non-Ohmic behavior to be 
associated with 'hidden' structural defects. The thickest 
samples demonstrated relatively sharp R(T) supercond- 
cuting transitions with the shape which can be qualita- 
tively understood by the model of thermally activated 
phase slips (TAPS)i^ and the inevitable mild inhomo- 
geneity of the samples 6 . However for the nanowires with 
diameters < 50 nm the width of the R(T) transition 
broadens well above the limits which can be explained by 
the TAPS model with a realistic set of parameters. For 
the thinnest samples with diameters < 30 nm the temper- 
ature dependence of the experimentally measured resis- 
tance R{T) is very weak and does not extrapolate to zero 
at T — > 0. By the standards of modern nanotechnol- 
ogy and microscopic analysis the structures are so 'large' 
that speculations about trivial inhomogeneity overlooked 
in SET, TEM and SPM microscopes could be ruled out 
with a high level of confidence. We associate the pro- 
nounced broadening of the R(T) transitions with quan- 
tum fluctuations of the order parameter - the so-called 
quantum phase slips (QPS^rij^. Additionally to the 
importance for the basic knowledge about nanoscale su- 
percondcutivity, the subject of quantum fluctuations is 
expected to lead to a new class of devices: quantum stan- 
dard of electric current^ 6 -, qubilj^l and various QPS-based 
system s 38 , 39 . 
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